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The correlation function of two identical particles - pions or kaons - interacting via Coulomb 
potential is computed. The particles are emitted from an anisotropic particle's source of finite 
lifetime. In the case of pions, the effect of halo is taken into account as an additional particle's 
source of large spatial extension. The relativistic effects are discussed in detail. The Bowler- 
Sinyukov procedure to remove the Coulomb interaction is carefully tested. In the absence of halo 
the procedure is shown to work very well even for an extremely anisotropic source. When the halo 
is taken into account the free correlation function, which is extracted by means of the Bowler- 
Sinyukov procedure, is distorted at small relative momenta but the source parameters are still 
Oh' correctly reproduced. 
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I. INTRODUCTION 

The correlation functions of two particles with small relative momenta provide information about space-time char- 
acteristics of particle's sources in high-energy nucleus-nucleus collisions, see the review articles [3, 0, Within the 
standard femtoscopy, one obtains parameters of a particle's source, comparing the experimental correlation functions 
to the theoretical ones which are calculated in a given model. Such an analysis can be performed for pairs of non- 
identical or identical particles. In the former case, the correlation appears due to inter-particle interaction while in the 
latter one the interaction is combined with the effects of quantum statistics. Since we usually deal with electrically 
charged particles, observed two-particle correlations are strongly influenced by the Coulomb interaction. The effect of 
the Coulomb force is usually eliminated from experimental data by means of the so-called Bowler-Sinyukov procedure 
[B|, 0] . And then, the correlation function, which is obtained in such a way from experimental data, is compared to the 
theoretical correlation function of two non-interacting particles. The comparison provides parameters of the source 
of particles. 

\ The femtoscopy was applied to a large volume of experimental data on nucleus- nucleus collisions at SPS energy as 
summarized in [7j. The spatial size of particle's sources appeared to be comparable to the expected size a fireball 
Q\ | created in nucleus-nucleus collisions while the emission time of particles was significantly shorter. It was predicted 
that at RHIC energies the emission time would be significantly longer due to the long lasting hydrodynamic evolution 
of the system created at the early stage of nucleus- nucleus collisions [1, Q . To a big surprise the experimental data 
obtained at RHIC [l(| HH, EH EH show a very little change of the space-time characteristics of a fireball when compared 
to the SPS data. In particular, the emission time of particles appeared to be as short as 1 fm/c. Because of this 
surprising result, which is now known as the 'HBT Puzzle' [Til. Ila]. a reliability of the femtoscopy was questioned. 

As an alternative to the standard femtoscopy, the method of imaging E3] was developed. Within this method 
one obtains the source function not referring to its specific parametrization but directly inverting the correlation 
function. The procedure of inversion takes into account the effect of quantum statistics as well as that of inter- 
particle interaction. The one-dimensional and three-dimensional imaging was successfully applied to experimental 
data, see [HI, EH and [13, HH, HH, respectively. The method provides essentially model independent information 
on the source space-time sizes but modeling is still needed to deduce the emission time which is coupled to spatial 
parameters of the source. Therefore, the imaging has not much helped to resolve the 'HBT puzzle'. 

Very recently it has been shown that hydrodynamic calculations can be modified to give quite short emission times 
of produced particles [23[ , see also 0, 25] . Specifically, the initial condition needs to be changed to speed up formation 
of the transverse collective flow and the first order phase transition from quark-gluon plasma to hadron gas should 
be replaced by the smooth cross-over. Another solution of the 'HBT Puzzle' assumes an incomplete equilibrium of 
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quark-gluon plasma [26[. Although the 'HBT Puzzle' is resolvable now, if not resolved in [H, [H, [H, H^l i it is still of 
interest to quantitatively check the femtoscopy method, to be sure that experimentally obtained source parameters 
are indeed reliable. This is the aim of our study which is mainly focused on the Coulomb effects. Our preliminary 
results were presented in [27j but, unfortunately some errors appeared in this publication. 

The Bowler-Sinyukov correction procedure, which is used to eliminate the Coulomb interaction from the experimen- 
tal data, assumes that the Coulomb effects can be factorized out. The correction's factor is calculated for a particle's 
source which is spherically symmetric and has zero lifetime. We examine the procedure applying it to the computed 
Coulomb correlation functions of identical particles coming from aniso trop ic sources of finite lifetime. Azimuthally 
asymmetric sources, which appear in azimuthally sensitive femtoscopy USUI]' are also studied. We treat the com- 
puted Coulomb correlation functions as experimentalists treat the measured correlation functions. Thus, we extract 
the correlation function which is supposed to be free of the Coulomb interaction. However, in contrast to the situation 
of experimentalists we know actual parameters of particle sources which can be compared to the extracted ones. Our 
analysis is somewhat similar to that presented in [30| where the Coulomb correlation functions were computed by 
means of Monte Carlo event generator and it was claimed that the procedure of removal of the Coulomb effects works 
well. Our analysis is more detailed and it is based on mostly analytical calculations. 

The correlation function of two identical non-interacting bosons is expected to be equal to 2 for vanishing relative 
momentum of the two particles. The correlation functions extracted from experimental data by means of the proce- 
dure, which is supposed to remove the Coulomb interaction, do not posses this property. The correlation function 
at zero relative momentum is significantly smaller than 2. This fact is usually explained referring to the concept 
of halo [31| . It assumes that only a fraction of observed particles comes from the fireball while the rest originates 
from the long living resonances. Then, we have two sources of particles: the fireball and the halo with the radius 
given by the distance traveled by long living resonances. The complete correlation function, which includes particles 
from the fireball and the halo, equals 2 at exactly vanishing relative momentum. However, the correlation of two 
particles coming from the halo occurs at a relative momentum which is as small as the inverse radius of the halo. 
Since experimental momentum resolution is usually much poorer and such small relative momenta are not accessible, 
the correlation function is claimed to be less than 2 for effectively vanishing relative momentum. We carefully study 
the effect of halo and, in particular, we test how the Bowler-Sinyukov correction procedure works in the presence of 
halo. 

We discuss in detail how to compute the Coulomb correlation functions. We pay particular attention to relativistic 
effects which, in our opinion, are not clearly exposed in literature. We start with the nonrelativistic Koonin formula [32[ 
because of its rather transparent physical meaning. A more formal derivation of the correlation function, which follows 
the studies [H, [34], [Hj], is sketched in the Appendix A. The Koonin formula expresses the correlation function through 
the nonrelativistic wave function of two particles of interest. Since the observed correlation functions are significantly 
different from unity only for small relative momenta when the relative motion of particles is nonrelativistic, it is 
legitimate to use the nonrelativistic wave function in the center-of-mass frame of two particles. However, it requires an 
explicit transformation of the source function to the center-of-mass frame. It should be mentioned that transformation 
properties of nonrelativistic wave function under a Lorentz boost are not well understood. Only recently it has been 
shown using the Bethe-Salpeter equation that the hydrogen atom wave function experiences the Lorentz contraction 
[36| under the Lorentz boost. Therefore, we perform the calculation in the center-of-mass frame of the pair and then 
we transform the correlation function to the source rest frame. 

Throughout our whole analysis the source function is of the Gaussian form. Such a choice has several advantages. 
First of all, the Gaussian source functions are often used to describe experimental data. Actually, the i maging method 
QIKItJ shows that non-Gaussian contributions to the source functions are at a percent level [H, QjJ [2(| HJ H2] • There 
are also pure theoretical advantages of the Gaussian parameterization. When the single particle source function is 
Gaussian, the so-called relative source function is Gaussian as well. Since the Gaussian source function has a simple 
Lorentz covariant form, Lorentz transformations can be easily performed. Due to the two features of the Gaussian 
source functions, our calculations are mainly analytical which in turn allowed us, in particular, to carefully study 
relativistic effects mentioned above. 

The Gaussian parameterization we use has an important disadvantage - the fireball expansion is entirely neglected. 
The study of source expansion, however, goes beyond the scope of our analysis. We address in this paper a specific 
question whether the Bowler-Sinyukov procedure properly removes Coulomb effects from the correlation functions. For 
this purpose we compute the Coulomb correlation function with the Gaussian source, we apply the Bowler-Sinyukov 
procedure and we check how accurately the free correlation function, which is also computed with the Gaussian 
source, is reproduced. It is certainly of interest to study how the fireball's expansion influences the Coulomb effects. 
Before that, however, one should systematically analyze to what extend the expanding fireball can be represented by 
a Gaussian source. For this reason we do not discuss the interplay of Coulomb interaction and fireball expansion. 
Actually, the problem cannot be addressed using the computational methods we developed. 

Throughout the paper we use natural units, where c — % — 1, and our metric convention is (+, — , — , — ). 
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FIG. 1: The -kit Coulomb correlation functions C(q , 0,0) (left panel), C(0,g s ,0) (central panel) and C(0, 0, q{) (right panel) 
as functions of q , q s or qi, respectively, for three values of the emission time r = 1,2,3 fm. The remaining parameters are: 
R x = 4 fm, R y = 4 fm, R z = 6 fm, and v = (0.8, 0, 0). 



II. DEFINITION 



The correlation function C(pi, P2) of two particles with momenta pi and P2 is defined as 

_ dN 

C(pi,p 2 ) 



dp 1 dp 2 



dN dN 
dpi dp 2 



(1) 



where dp d ^ P2 and ^~ is, respectively, the two- and one-particle momentum distribution. The correlation function 
can be written down in a Lorentz covariant form 



C(pi,p 2 ) 
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(2) 



' dp 2 



where E4¥- is the Lorentz invariant distribution. 

d°p 



The covariant form ([2]) shows that the correlation function is a Lorentz scalar field which can be easily transformed 
from one reference frame to another. If the particle four-momenta, which are on mass-shell, transform as Pi — > p'j 
with i = 1,2, the transformed correlation function equals 

C'(pi(Pi),P 2 (p 2 ))-C < (p 1 ,p 2 ). 



III. NONRELATIVISTIC KOONIN FORMULA 

Within the Koonin model [HI], the correlation function C can be expressed in the source rest frame as 

C(pi,p 2 ) - J d 3 r 1 dt 1 d 3 r2dt 2 D(t 1 ,r 1 ) D(t 2 ,r 2 ) |*(ri, r' 2 )| 2 , (3) 

where v\ = ri+vrfi, r' 2 ) is the wave function of the two particles and D(t, r) is the single-particle source function 

which gives the probability to emit the particle from the space-time point (t, r). The source function is normalized as 

J d 3 rdtD(t,r) = 1 . (4) 

After changing the variables r' «-> r, the correlation function can be written in the form 

C(pi, p 2 ) = J d 3 r 1 dt 1 d 3 r 2 dt 2 D(t 1 ,r 1 - viii) D(t 2 , r 2 - v 2 t 2 )|*(n, r 2 )| 2 . 

Now, we introduce the center-of-mass coordinates 

r = r 2 -ri, R = -^(miri + m 2 r 2 ), 

t = t 2 -h, T = jj{m 1 t 1 +m 2 t 2 ), 

q= Ti(m 2 pi -mip 2 ), P = pi+p 2 , 
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FIG. 2: The KK Coulomb correlation functions C(q o ,0, 0) (left panel), C(0,q 3 ,0) (central panel) and C(0,0, qi) (right panel) 
as functions of q , q s or qi, respectively, for three values of the emission time r = 1,2,3 fm. The remaining parameters are: 
R x = 4 fm, R y = 4 fm, R z = 6 fm, and v = (0.8, 0, 0). 



where M = mi + m 2 . Using the center-of-mass variables, one gets 



C(q)= / d 3 rZ> r (r)|^ q (r)| 2 



(5) 



where the 'effective relative' source function is defined as 



D r (r) ee J dtD r (t,r-vt) , (6) 
and the 'relative' source function is expressed through the single-particle source function in the following way 

A.(i,r) = Jd 3 RdTD(T~^R-^r)D(T+^t,R+^r). (7) 
We note that due to the normalization the functions D r (r) and D r (r — vt,t) are also normalized 

J d 3 rD r {r) = J d 3 rdt D r (t,r) = 1 . (8) 



To get Eq. ([5]), the wave function was factorized as 

*0i,r 2 ) 



iPR 



^q(r) 



with y q (r) being the wave function of the relative motion in the center-of-mass frame. Deriving Eq. ([5]), it has been 
assumed that the particle velocity, which enters the effective source function, is the same for both particles. Thus, we 
have assumed that vi = v 2 = v which requires, strictly speaking, q = 0. However, one observes that |vi — v 2 | <C |vj| 
if |q| <C n\pi\/mi where fi = mim^/M. Thus, the approximation vi w v 2 holds for sufficiently small particle's 
momenta in the center-of-mass frame. It should be stressed that the dependence of the correlation function on q is 
mostly controlled by the dependence of the wave function on q which is not influenced by the above approximation. 
We choose the Gaussian form of the single-particle source function D(t, r) 



D(t,r) 



1 



4tt 2 R x R v R z t 



exp 



2r 2 2i?2 2R 2 y 2RI 



(9) 



where r = (x, y, z) and the parameters t, R x , R y and R z characterize the lifetime and sizes of the source. Specifically, 
the parameters r and R x give, respectively, 

t 2 = (t 2 ) ee J d 3 rdtt 2 D(t,r) , R\ = (x 2 ) = J d 3 r dt x 2 D(t,r) . 



The relative source function computed from Eq. |(7J) with the single-particle source © is 



D r (t,r) 



1 



16it 2 R x R v R z t 



exp 



4t 2 4i?2 4i?2 AR 2 



(10) 
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FIG. 3: The mr Coulomb correlation functions C(g o ,0, 0) (left panel), C(0, q 3 ,0) (central panel) and C(0,0,qi) (right panel) as 
functions of q , q s or qi, respectively, for various halo contributions. The halo is spherical with Rh == 40 fm and th = 0. The 
fireball parameters are R x = 4 fm, R y — 4 fm, R z — 6 fm, r = 3 fm. The pair velocity is v = (0.8, 0, 0). 



We note that the particle's masses, which are present in the definition {?]), disappear completely in the formula (| 10[) . 
This is the feature of the Gaussian parameterization ((9]). 

In the case of non-interacting identical bosons, the two-particle symmetrized wave function is 

tf(ri,T 2 ) = -L[ e l Pi r i+P^2 + e ip 2 rx+pir 2 j _ [ e iqr + e -iqr] e iPR_ 

v 2 v 2 

It gives the modulus square of the wave function of relative motion |y> q (r)| 2 = 1 + cos (2qr) which in turn provides 
the correlation function equal to 

C(q) = 1 +exp [- 4(r 2 (qv) 2 + R 2 x q 2 x + R 2 y q 2 y + R 2 z q 2 z )] , (11) 

where q = (q x ,q y ,q z ). We note that the 'cross terms' such as q x q z , which are discussed in [U, do not show up, 
as the source function ([9]) obeys the mirror symmetry D(t,r) = D(— £, — r). We also note that q often denotes the 
relative momentum pi — p 2 not the momentum in the center-of-mass frame, which for equal mass nonrelativistic 
particles equals ^(pi — P2), and then, the factor 4 does not show up in the correlation function (fTT|) of identical free 
bosons. However, we believe that using the momentum in the center-of-mass frame is physically better motivated, as 
the center-of-mass variables naturally appear when the center-of-mass motion is separated from the relative one. 



IV. RELATIVISTIC FORMULATIONS 



There are two natural ways to 'relativize' the Koonin formula ([3]). The first one provides an explicitly Lorentz 
covariant correlation function but it is applicable only for the non-interacting particles. The second one holds only in 
a specific reference frame but it is applicable for interacting particles as well. Below, we consider the two methods. 
We start, however, with the discussion of the Lorentz covariant form of the source function. 



A. Lorentz covariant source function 



Because of its probabilistic interpretation, the source function transforms under Lorentz transformation as a scalar 
field. Therefore, the covariant form of the Gaussian parameterization of the source function © is written as [38| 



D(x) - exph-^A^], (12) 

where is the position four-vector and is the Lorentz tensor characterizing the source which in the source rest 
frame is 
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FIG. 4: The correction factor K(q) divided by the Gamov factor G(q) as function of q for source radii R = 1, 3, 5 fm. The left 
panel is for pions and the right one for kaons. The solid lines and squares represent, respectively, the exact formula (|35|l and 
the approximate one (|36l) . 



The source function as written in Eq. (|12[) obeys the normalization condition not only for the diagonal matrix A 
but for non-diagonal as well. 

The source function (|12p is evidently the Lorentz scalar that is 



V dctA' r 1 , ,,,,,, ,. V detA . 1 , ,,,, . „, . 
D W) = ex P i~2< A M <\ = ~^2~ ex P [-2^ A Xv] = D{x) 



At: 2 



(14) 



where x'^ = L"x u and A' Miy = L fl a .A ap L p l/ with being the matrix of Lorentz transformation. We note that 
detA' = detL detA detL" 1 = detA. 

The covariant relative source function (|10p is given by 



VdetA 1 
A-(» = — „ exp ~7^ A ^ • 



(15) 



B. Explicitly covariant 'relativization' 

As follows from Eq. ((2]), the correlation function is a Lorentz scalar. Therefore, the Koonin formula (J3j) can be 
'relativized' demanding its Lorentz covariance. Let us write the formula as 

Cfa,P2)= f d 4 x 1 ct k x 2 D(x 1 )D(x 2 )\^{x 1 ,x 2 )\ 2 , (16) 

where pi and Xi is, respectively, the four-momentum and four-position. Since the source function D(x) and the four- 
volume element d 4 Xi are both the Lorentz scalars, the whole formula (|16p is covariant if the wave function ^>(xi, x 2 ) 
is covariant as well. In the case of non- interacting bosons the relativistic wave function ^(xi, x 2 ) is 

v2 

As the function depends on the scalar products of two four-vectors, it is the Lorentz scalar. We note that the function 
(|17p depends on two time arguments. 

Our further considerations are limited to pairs of identical particles and thus, we introduce the relative coordinates: 

x = x 2 — xi, X = \(xi +x 2 ), 

(18) 

q= 2 (pi-P2), P = pi+p 2 - 

In this section and in Appendix A q is the four-vector (qo, q) but in the remaining sections q = |q|. Hopefully, it will 
not cause any confusion. 

We note that in the nonrelativistic treatment the three- vectors r and q, which are given by the four-vectors x = (t, r) 
and q — (go, q), correspond to the inter-particle separation and the particle's momentum in the center-of-mass of the 
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pair. This is, however, not the case in the relativistic domain. To get the center-of-mass variables, the four-vectors 
need to be Lorentz transformed. We also note that qo = qv which is proven as .45] 
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2 yj m 2 + p 2 + yjrn 2 
With the variables (TS)) . the wave function (fTT|) equals 

9{x,X) = -^(e iqx + e- 



' a/to 2 + pf + sjm 2 + 



p! 



qv 



PX 



(19) 



and the correlation function is found in the form 

C(g) = l + exp[-%(A^)-V] , 

which is explicitly Lorentz covariant. For the source matrix (|13p . the correlation function equals 

C(q) = 1 + exp [ - A{q 2 T 2 + q 2 x R 2 x + q 2 y R 2 y + q 2 z R 2 z )] . 

Since qo = qv, the correlation function (flT))) exactly coincides with the nonrelativistic expression (fTTj) . This coincidence 
is not completely obvious as the time variables enter differently in the Koonin formula ([3]) and in the covariant one 

m. 

Let us consider the correlation function in the center-of-mass frame of the particle pair. We assume that the 
velocity of the center-of-mass frame in the source rest frame is along the axis x. Then, v = (v,0, 0) and qo = q x v. 
The correlation function (| 19[) . which holds in the source rest frame, equals 

1 + exp [ - A((v 2 t 2 + Rl)ql + R 2 y q 2 y + R 2 z q 2 z )} . 



C(q) 



(20) 



As seen, the effective source radius in the direction x is y/ R 2 + v 2 t 2 . We now transform the source function to the 
center-of-mass frame where the quantities are labeled with the index *. The center-of-mass source matrix (|13[) . which 
is computed as 
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FIG. 6: The KK free correlation function C(q , 0,0) (left panel), C(0,q s ,0) (central panel) and C(0, 0, qi) (right panel) as 
a function of q , q s or qi, respectively. The fireball parameters are R x = 4 fm, R y = 4 fm, R z — 6 fm, r = 3 fm and the 
pair velocity is v = (0.8,0,0). The extracted 'free' functions are represented by the squares and the expected free correlation 
functions by the solid lines. 



Then, the correlation function in the center-of-mass frame is found to be 

C7(g.) = 1 + expMg^ArrV,] = 1+ exp [ - 4( 7 2 ( V 2 r 2 + R 2 x )ql x + R 2 y ql y + R 2 z ql z )] . 



(22) 



As seen, the effective source radius along the direction of the velocity is elongated, not contracted as one can naively 
expect, by the factor 7. 



C. Non-covariant relativization 



The quantum mechanical description of two relativistic interacting particles faces serious difficulties. The problem 
is greatly simplified when the relative motion of two particles is nonrelativistic (with the center-of-mass motion 
being fully relativistic). Since the correlation functions usually differ from unity only for small relative momenta of 
particles, it is reasonable to assume that the relative motion is nonrelativistic. We further discuss the correlation 
functions taking into account the relativistic effects of motion of particles with respect to the source but the particle's 
relative motion is treated nonrelativistically. In such a case, the wave function of relative motion is a solution of the 
nonrelativistic Schrodinger equation. Thus, we compute the correlation function directly from the Koonin formula ([3]) 
but the computation is performed in the center-of-mass frame of the pair. For this reason we first transform the source 
function to this frame and then, after performing the integrations over x\ and X2, we transform the whole correlation 
function to the source rest frame. We stress here that according to the definition @ the correlation function is the 
Lorentz scalar. 

As already noted, we compute the correlation function in the center-of-mass frame of the pair using the relative 
variables (fT8|) . The correlation function thus equals 



C(q*) = J d 3 r*dt* D r (i*,r*) |v? q „ (r»)| 2 , 



(23) 



where -D r (i», r») is the relative source function (j 1 5[) and ipq t (r„) is, as previously, the nonrelativistic wave function of 
relative motion. The note here that v* = by definition. The formula (|23[) can be rewritten as 



C(q*) = J d 3 ^ D r (r*) \(PdA r *)\' 



where 



D r (r») = / dt*D r (t*,r*) = — 
J 8tt 



3 / 2 y/-f 2 (R 2 x +V 2 T 2 )R y R, 



■ exp 



J 2 (Rx + V 2 T 2 ) 



Rl 



4.) 



(24) 



(25) 



for v = (v, 0, 0). One easily checks that the free correlation function, which follows from Eq. ([23]) or Eq. ([24]) . exactly 
coincides with the formula (|22|) . To get the correlation function in the source rest frame, one performs the Lorentz 
transformation (the correlation function as defined by Eq. (J5|) is a Lorentz scalar) and obtains the formula ([2H|) . Thus, 
the two ways of 'relativization' give the same result for non-interacting particles. This is not quite trivial as the time 
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FIG. 7: The -kit Coulomb correlation functions C(q , 0,0) (left panel), C(0,g s ,0) (central panel) and C(0, 0, qi) (right panel) 
as functions of q , q 3 or qi, respectively, for the azimuthally asymmetric source. The parameters are: R x — 2 fm, R y = 4 fm, 
R z = 6 fm, t = 1 fm and the pair velocity is v = (0.8, 0, 0). The extracted 'free' functions are represented by the squares while 
the expected free correlation functions correspond to the solid lines. 



dependence of the Koonin formula (f2"3"| and of the explicitly covariant one fTO)) look rather different. Unfortunately, 
we do not know whether the equivalence of the two 'relativization' schemes holds for interacting particles, as the 
covariant 'relativization' is known only for free particles. 



V. COULOMB CORRELATION FUNCTIONS 



In this section we compute, using Eq. (f2"4")l . the correlation functions of pairs of identical pions or kaons interacting 
due to the Coulomb force. The calculations arc performed for the anisotropic Gaussian source of finite emission time 
(fT2l [TU]) . We use the Bertsch-Pratt coordinates [4j], |42j out, side, long. These are the Cartesian coordinates, where 
the direction long is chosen along the beam axis (z), the out is parallel to the component of the pair momentum P 
which is transverse to the beam. The last direction - side - is along the vector product of the out and long versors. 
So, the vector q is decomposed into the q , q s , and qi components. If the particle's velocity is chosen along the axis 
x, the out direction coincides with the direction x, the side direction with y and the long direction with z. We note 
that the correlation function of two identical free bosons in the Bertsch-Pratt coordinates in the source rest frame is 

<7(q) = 1 + exp [ - 4(q 2 R 2 + q 2 s R 2 s + q 2 R 2 

where R a = \J R 2 . + v 2 t 2 , R s = R y and Ri = R z . As seen, the source lifetime r is mixed up with the size parameter 
R x . Although experimentalists usually use the parameters R Q , R s , Ri, we use them together with R x , R y , R z and r, 
as the lifetime r naturally enters theoretical formulas. Since the velocity of the pair is chosen along the axis x, we 
always have R s = R y and Ri — R z . 

The effect of Coulomb interaction in femtoscopy can be treated analytically or almost analytically under some 
simplifying approximations 43]. However, we are interested in the exact Coulomb correlation functions. Therefore, 
we use the exact wave function. In the case of two non-identical particles interacting due to repulsive Coulomb force, 
the nonrelativistic wave function is well known to be [39j 

p q (r) - e~% r(l + d) e^ ^(-i^l, i{qr - qr)) , (26) 
q q 

where q = |q| and rj^ 1 is the Bohr radius which for pairs of pions and kaons equals rj^ 1 = 388 fm and rf^ = 110 
fm, respectively; %Fi denotes the hypergeometric confluent function. When one deals with identical bosons, the wave 
function </3 q (r) should be symmetrized and the modulus of the symmetrized Coulomb wave function equals 



1 



■ 7 1 1 -I _„\M2 i I 777 / -V 



Pq(r)| J = ^G(q) | 1 F 1 H^,l,z( <Z r~qr))| i + | 1 ^ 1 (- l ^,l,z( (Z r + qr))| 2 (27) 



q q 



2Re(e 2 "i r iFi(-*2 1, i(qr - qr)) if?H~, 1. Kv + w)) 



where G(q) is the so-called Gamov factor defined as 

c 



G{q) = ^ - r L 7 - (28) 
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FIG. 8: The irn Coulomb correlation functions C(q , 0,0) (left panel), C(0,g s ,0) (central panel) and C(0, 0, qi) (right panel) 
as functions of q , q 3 or qi, respectively, for the azimuthally asymmetric source. The parameters are: R x — 4 fm, R y — 2 fm, 
R z = 6 fm, t = 1 fm and the pair velocity is v = (0.8, 0, 0). The extracted 'free' functions are represented by the squares while 
the expected free correlation functions correspond to the solid lines. 



One may wonder whether a multi-particle environment, which occurs in the final state of relativistic heavy- ion 
collisions, influences the Coulomb potential of the two particles of interest. It should be remembered, however, that 
the particles are correlated at small relative momenta and thus, they fly with similar velocities. Consequently, after 
the time comparable to the source size, the particles with small relative velocity appear to be effectively isolated from 
the rest of many-particle system. Therefore, the effect of screening of Coulomb potential is expected to be negligible. 
This qualitative argument is confirmed by the calculations presented in [40j. 

Substituting the modulus l|27p and the source function (|25p into Eq. (|2~4")) , one finds the correlation function in 
the center-of-mass frame which is further transformed to the source rest frame. In Figs. [1] and [2] we show the 
correlation functions C(g o ,0,0), C(0,q s ,0) and C(0, 0, qi) of identical pions and kaons, respectively. The calculations 
are performed for the following values of the source parameters: R x — 4 fm, R y — 4 fm, R z — 6 fm, and r = 1, 2, 3 fm. 
The velocity of the particle's pair with respect to the source equals v = 0.8 and it is along the axis x. As seen, the 
correlation functions of pions and kaons differ sizably due to the different Bohr radii of the two systems. The most 
visible difference appears for the function C(q o ,0, 0). 



VI. THE HALO 

As mentioned in the introduction, the halo [31( was introduced to explain the fact that, after removing the Coulomb 
effect, the experimentally measured correlation functions are smaller than 2 at vanishing relative momentum. The 
idea of halo assumes that only a fraction / (0 < / < 1) of particles contributing to the correlation function comes 
from the fireball or core while the remaining fraction (1 — /) originates from long living resonances. Then, we have 
two sources of the particles: the small one - the fireball or core - and the big one corresponding to the long living 
resonances. The single-particle source function has two contributions 

D(t, r) = / D f (t, r) + (1 - /) D h (t, r) , (29) 

where Df(t,r) and Dh(t,r) represent the fireball and halo, respectively. For non-interacting identical bosons, the 
correlation function is 

C(q) - 1 + fe- iR W + (1 - f) 2 e~ iR2 ^ 2 + 2/(1 - f) e ^{R}+KW , ( 30 ) 

where both the fireball and halo are assumed to be spherically symmetric sources of zero lifetimes; Rf and Rh are 
the radii of, respectively, the fireball and the halo. If Rh is so large that R^ 1 is below an experimental resolution of 
the relative momentum q, the third and fourth term of the correlation function (|30[) are effectively not seen, and one 
claims that C(q = 0) = 1 + A with A = / 2 < 1. 

We have included the halo in our calculations of the 7T7t Coulomb correlation functions. Since the halo represents 
pions from resonances, the source function of the halo, which was carefully modeled in |44j], is approximately of the 
exponential form. In our calculations, however, the halo source function, as other source functions we use, is of 
the Gaussian form for the reasons explained in the introduction. Our simplified treatment of the halo seems to be 
harmless, as the halo influences the correlation function only for q of the order R^ 1 which are experimentally hardly 
accessible. 



11 




q [MeV/c] q s [MeV/c] q| [MeV/c] 



FIG. 9: The mr and KK Coulomb correlation functions C(q , 0,0) (left panel), C(0, g s ,0) (central panel) and C(0,0,<#) (right 
panel) as functions of q , q s or g;, respectively, for the extremely anisotropic source. The parameters are R x — R y — 0, 
R z = 6 fm, t = 3 fm and the pair velocity is v = (0.8, 0, 0). 



Our exemplary results are shown in Figs. [3] for several values of A. The fireball is anisotropic with R x = 4 fm, 
R y = 4 fm, R z = 6 fm, r = 3 fm and v — 0.8; the halo is spherically symmetric, its radius is Rh = 40 fm (as suggested 
in [3l|) and its lifetime vanishes. 

In principle, a finite spatial extension of the halo implies a finite duration of pion emission. However, when the 
finite lifetime of the halo is taken into account, the size of the halo in the out direction increases, and the correlation 
function observed in this direction is influenced at even smaller momenta than those in the side and long directions. 
In other words, neglecting the finite lifetime of the halo, its effect on the correlation function in the out direction is 
overestimated not underestimated. We return to this point at the end of Sec. [X] 



VII. COULOMB CORRECTION WITHOUT HALO 



As mentioned in the Introduction, the Coulomb effect is usually subtracted from the experimentally measured 
correlation functions by means of the Bowler-Sinyukov procedure. We first note that the Coulomb effect is far not 
small and thus the method to subtract the Coulomb effect should be carefully tested. 

In the absence of halo the Bowler-Sinyukov procedure assumes that the Coulomb effect can be factorized out, that 
is the correlation function can be expressed as 

C(q) = K{q) C free (q) , (31) 

where Cf rec (q) is the free correlation function and K(q) is the correction factor which can be treated as the Coulomb 
correlation function of two nonidentical particles of equal masses and charges. The function is, however, rather 
unphysical as the pair velocity vanishes even so the calculation is performed in the rest frame of the source where the 
source is assumed to be symmetric and of zero lifetime. The correction factor K(q), which is described in detail in 
the Appendix to the paper [30l| . is computed as 

K(q) = G(q) /dVA-MliF^-^M^r-qr))! 2 , (32) 
J 1 

where G(q) is the Gamov factor ([28]) and D r (r) describes the spherically symmetric Gaussian source of zero lifetime 
and of the 'effective' radius R — y/ {R 2 , + R% + Ri)/3 where R — \JR% + v 2 t 2 , R s — R y and Ri — R z are the 
femtoscopic radii obtained from the extracted free correlation function. Experimentally R , R s and Ri are found 
fitting the measured correlation function C(q) with K(q) Cf roc (q). In our theoretical analysis, R x , R y , R z and r are 
the actual source parameters which enter the source function (J9j> . 

Using the parabolic coordinates £ + = r + z, £_ = r — z and the azimuthal angle <f), the relative source function of 
isotropic Gaussian source of zero lifetime is 

^+,^) = ^exp(-%^), (33) 
which substituted into Eq. (|32p gives 



iFl(--,l,i5f-) 
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FIG. 10: The nn and KK 'free' correlation functions C(q , 0,0) (left panel), C(0,q s ,0) (central panel) and C(0,0, qi) (right 
panel) as functions of q , q s or qi, respectively, for the extremely anisotropic source. The parameters are: R x = R y = 0, 
R z = 6 fm, t = 3 fm and the pair velocity is v = (0.8, 0, 0). The extracted 'free' functions are represented by the squares for 
pions and by triangles for kaons; the expected free correlation functions correspond to the solid lines. 



The trivial integral over <f> has been performed in Eq. (|34[) . Since the confluent hypergeometric function does not 
depend on £ + , the integral over £ + can be easily performed and one obtains 



K{q) = 



G(q) 
2tt 1 / 2 R 



q 



exp 



e 



16R 2 



(35) 



where the integral over £_ is usually computed numerically. However, observing that the source size is always much 
smaller than the Bohr radius of the particles of interest, one derives the approximate expression of the hyperbolic 
confluent function (|B9|) which is discussed in Appendix B. With the formula (|B10|) the integration can be performed 
analytically and the correction factor equals 



K(q) = G(q) 



1 



8r]R 



•Fn 



13 3 

2' '2' 2' 



-Aq 2 R 2 



(36) 



In Fig. 0] we show the correction factor K{q) for pions and kaons computed from the exact formula (j3"2"|) and the 
approximate one (f36|) . To make the difference more visible (note the vertical scale) the correction factor is divided 
by the Gamov factor which strongly varies with q. One sees that the approximation (|36p is very accurate for pions 
and it is less accurate for kaons. For this reason the expression (|3"6"|) is used only for pions. Fig. [3] also shows that the 
correction factor is heavily dominated by the Gamov factor that is K(q)/G(q) differs very little from unity. 

Once we are able to compute the exact Coulomb correlation functions for an anisotropic source of finite lifetime, 
we can test whether the free correlation function obtained by means of the Bowler-Sinyukov equation (|31[) properly 
reproduces the actual free correlation function. 

The free correlation functions, which are obtained using Eq. (|3"Tj) , are shown in Fig. [5] for pions in Fig. [5] for kaons. 
The extracted functions are compared to the expected correlation functions of noninteracting bosons for the given 
source. As seen, the free correlation function is almost exactly reproduced in the case of pions while in the case of 
kaons the reproduction is less accurate. Similar results are found as long as the source radii are much smaller than 
the Bohr radius of particles of interest. 



VIII. AZIMUTHALLY SENSITIVE FEMTOSCOPY 



In the previous sections we discussed the particle sources of cylindrical (azimuthal) symmetry (R x = R y ). The 
sources created in non-central collisions are not azimuthally symmetric but the symmetry is usually restored due to 
the averaging over impact parameter orientation. The cylindrically asymmetric sources are observable, if the reaction 
plane is determined. The azimuthally sensitive femtoscopy was developed [28l . [2^ | and when applied to experimental 
data it showed an expected dependence of the source radii on the emission angle with respect to the reaction plane. 
Since the Coulomb effects were removed from the data by means of the Bowler-Sinyukov procedure in the experimental 
studies [H,[2§], we test in this section the procedure for the case azimuthally asymmetric sources. 

We do not study a full dependence of the correlation function on the azimuthal emission angle but we consider two 
extreme cases. We stick to our convention that particles are always emitted along the axis x but R x ^ R y . In Fig. [7] 
we show the pion correlation functions C(g o ,0,0), C(0,q s ,0) and C(0, 0,qi) for R x = 2 fm, R y = 4 fm and in Fig. [5] 
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FIG. 11: The irn Coulomb correlation functions C(q o ,0, 0) (left panel), C(0, q 3 ,0) (central panel) and C(0, 0,qi) (right panel) 
as functions of q , q s or qi, respectively. The fireball and halo parameters are R x = 4 fm, R y — 4 fm, R z — 6 fm, r = 3 fm, 
Rh — 40 fm and the pair velocity is v = (0.8,0,0). The exact Coulomb correlation functions are shown with the squares and 
the Coulomb correlation functions fitted with the dilution formula (|4ip are represented by the solid lines. 

the case R x = 4 fm, R y = 2 fm is illustrated. The remaining parameters are R z = 6 fm, r = 1 fm, v = (0.8, 0, 0). As 
seen, the Bowler-Sinyukov procedure works very well in both cases. 



IX. EXTREMELY ANISOTROPIC SOURCE 



To establish limitations of the Bowler-Sinyukov procedure in the absence of halo, we have considered an extremely 
anisotropic source where R z is much larger than R x and R v as well as R z 3> r. The source function is found from 
Eq. (|25|) by taking the limits R y — ► 0, R x — > and r — > 0. Thus, one finds 



A-(r*) = 



1 



2n 1 / 2 R. 



■ exp 



ml 



(37) 



A paradoxical feature of this source function is that the information about the velocity v of the pair's center-of-mass 
frame with respect to the source has disappeared. Thus, we have the same source function in the source rest frame 
and in the center-of-mass frame of the pair. However, when we transform the correlation function from the pair 
center-of-mass to the source rest frame the pair velocity enters. The advantage of Eq. (|37|) is that the calculations 
can be performed almost analytically. 

Substituting the source function (|37[) into Eq. ([24]) . one finds 



ck,o,o) 
c(o, q ;,o) 



2G(q x ) f . . . rj . |2 --f^2 

rj} \ dz\ 1 F 1 {-i — ,l,iq x z)\ e 4R * , 
V™* Jo Qi 



2G( g ; 



G(q. 



z JO % 



2y/irR z J q z 



(38) 
(39) 
(40) 



-2 / dz 
lo 



Re[e- 2i "> 1 F 1 (-i^,l,2iq*z) ) I 



For pions the integration over z is performed using the approximate expression of the hypergeometric confluent 
function (IB9I). Thus, we find 



C(£,0,0) - 2G(q* x ) 
C(0,q;,0) = 2G{q* y ) 



^ 2 F 2 (i li^jf) 
V^F 2r2 4' i ' 2 ' 2' ^ zJ 

Jtx 2 2l 2' ' 2' 2' ^ zJ 
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FIG. 12: The irn Coulomb correlation functions C(q o ,0, 0) (left panel), C(0, q s ,0) (central panel) and C(0, 0,g;) (right panel) 
as functions of q , q s or qi, respectively. The fireball and halo parameters are R x = 4 fm, R y — 4 fm, R z — 6 fm, r = 3 fm, 
Rh — 40 fm and the pair velocity is v = (0.8,0,0). The exact Coulomb correlation functions are shown with the squares and 
the Coulomb correlation functions fitted with the Bowler-Sinykov formula (|42p are represented by the solid lines. 



, n (4q*R z f n n\ 2 F 2 (±,n + 1; §, §; -4q?R 2 z )\ 



" ^ ' ^ ^ ( 2 ") ! 2^(1,1; l 3 2 ;-MfRl) J ' 

The Bowler-Sinyukov procedure works very well for the correlation function computed with the extremely 
anisotropic source function (|3T[) with R z — 6 fm. To make the test of the procedure even more challenging, we 
have considered an extremely anisotropic source of finite lifetime. The Coulomb correlation functions of pions and 
kaons, which are computed for the source of R x — R y — 0, R z — 6 fm, r = 3 fm and v = 0.8, are shown in 
Fig. [5] The 'free' correlation functions, which are presented in Fig. 1101 are obtained from the correlation functions 
shown in Fig. [9] by dividing them by the correction factor K(q). The factor is computed for the averaged radius 
R = \J (Rq + Rf )/3- As seen, the Bowler-Sinyukov procedure works very well for both pions and kaons. In particular, 
C frcc (0, q 8 , 0) = C(0, q s ,0)/K(q) « 2 as expected. 



X. COULOMB CORRECTION WITH HALO 



The procedure to eliminate the Coulomb interaction is more complex when the halo is taken into account. We test 
two versions of the procedure which, following the STAR Collaboration [l3|, we call the 'dilution' method and the 
'proper Bowler-Sinyukov' one. The experimentally measured correlation functions C(q) are fitted as 



C(q) = (1 - A + XK(q)) [l + A(C frcc (q) - l) 
in the case of the dilution method and 



(41) 



C(q) = 1 — A + XK(q)C bee (q) , (42) 

in the case of the Bolwer-Sinykov method. 

The Coulomb correlation functions fitted according to the dilution (|4T|) and Bowler-Sinyukov (|4"2"|) formulas are 
shown in Fig. 1111 1121 respectively. The source parameters are given in the figures. The input parameters are those 
used in the computation of Coulomb correlation functions: the fireball parameters are R x = 4 fm, R y = 4 fm, 
R z — 6 fm and r = 3 fm; the halo is of zero lifetime of the radius Rh — 40 fm, and the pair velocity is v = (0.8, 0, 0). 
The output parameters are obtained from the fit. As seen, they only slightly deviate from the input parameters. Since 
the formulas (|41l I42p do not work at q < l/Rh, we perform the fit in the domain of q 0l q s or qi, respectively, bigger 
than 6 MeV. As seen the Coulomb correlation functions are fitted very accurately with both the dilution (|4"Tj) and 
Bowler-Sinyukov (|42[) formulas. 

The 'free' correlation functions extracted according to the dilution (|4"Tj) and Bowler-Sinyukov (T4"2"|) formulas are 
shown in Fig. 1131 1141 respectively. The expected free functions are also shown for comparison. It is important to 
note that the parameter A is assumed here to be known that is the actual value of A enters the formula (|4"1~T) or (142p . 
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FIG. 13: The tt-k free correlation functions C(q o ,0,0) (left panel), C(0,q s ,0) (central panel) and C(0,0, qi) (right panel) as 
functions of q , q s or qi, respectively. The fireball and halo parameters are R x = 4 fm, R y = 4 fm, R z = 6 fm, r = 3 fm, 
Rh = 40 fm and the pair velocity is v = (0.8,0,0). The 'free' correlation functions are extracted by means of the dilution 
method for several values of A. The actual free correlation functions are shown by the solid lines. 



As seen, the extracted correlation functions are distorted at small relative momenta and the distortions grow with A. 
However, the widths of the correlation functions are unaltered and so are the source parameters. 

The experimentally obtained 'free' functions, which are shown e.g. in Fig. 4 from (T3 |. do not reveal the dip at 
small q seen Fig. 13 and 14. We note, however, that the experimental correlations function in, say, out direction are 
not of the form C(q , 0, 0) but rather J dq s dqiC(q , q s ,qi) and the domain of truly small q is not seen. 

As mentioned at the end of Sec. IVI1 the model of halo should in principle include its finite lifetime. However, the 
finite lifetime increases the halo size in the out direction and the 'free' correlation function in out direction is distorted 
at even smaller momenta than the 'free' correlation functions in side and long directions. Therefore, our conclusions 
cannot be changed by taking into account a finite duration of pion emission from the halo. 



XI. CONCLUSIONS 



Let us summarize our study of the two-particle correlation functions. We have derived a relativistic generalization 
of the nonrelativistic Koonin formula. The calculations have been performed in the center-of-mass frame of the 
pair where a nonrelativistic wave function of the particle's relative motion is meaningful. It required an explicit 
transformation of the source function to the center-of-mass frame of the pair. Finally, the correlation function has 
been transformed to the source rest frame as a Lorentz scalar field. The Coulomb correlation functions of pairs of 
identical pions and kaons have been computed. The source has been anisotropic and of finite lifetime. For pions the 
effect of halo has been also taken into account. The source function has been always of the Gaussian form. 

Having the exact Coulomb correlation functions, the Bowler-Sinyukov procedure to remove Coulomb effect was 
tested. It was shown that the procedure works very well even for an extremely anisotropic source provided the halo 
is absent. For kaons small deviations are observed for a sufficiently large source. When the halo is included the pion 
correlation function are noticeably distorted for very small relative momenta but the source radii remain uninfluenced. 
Thus, we conclude that the Bowler-Sinyukov procedure, which at first glance does not look very reliable, appears to 
be surprisingly accurate. A possible interplay of Coulomb effects and fireball's expansion has not been studied here 
but our analysis shows that Coulomb effects are not sensitive to the source's shape as long as the characteristic source 
radius is much smaller than the Bohr radius of the particle's pair of interest. Then, the Bowler-Sinyukov procedure 
is expected to work well. 



APPENDIX A: DERIVATION OF CORRELATION FUNCTION 

We sketch here the derivation, which is discussed in detail in [H, H3, [35[ , of the correlation function of two identical 
interacting bosons. Under rather general conditions, the correlation function as defined by Eq. JT]) can be written 
down as 

CiPuPi) = J d 4 xid 4 x[d 4 X2d 4 x' 2 p(xi,X2;x[,X2)^p u p 2 (xi,x 2 )y* uP2 (x' 1 ,x' 2 ) , (Al) 

where p{x\, x 2 ; x\, x' 2 ) is the properly normalized coordinate space density matrix describing production process of 
the two particles and ^ Pl , P2 (xi, x 2 ) is the Bethe-Salpeter amplitude; x\, x%, x' 1i x' 2 , pi, p 2 are all four-vectors. 
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FIG. 14: The nn 'free' correlation functions C(q o ,0,0) (left panel), C(0,q a ,0) (central panel) and (7(0,0,5;) (right panel) as 
functions of q , q s or qi, respectively. The fireball and halo parameters are R x — 4 fm, R y = 4 fm, R z = 6 fm, r = 3 fm, 
Rh = 40 fm and the pair velocity is v = (0.8, 0, 0). The 'free' correlation functions are extracted by means of the Bowler-Sinykov 
method for several values of A. The actual free functions are shown by the solid lines. 



To separate the relative from center-of-mass motion, one uses the variables (118)) and expresses the Bethe-Salpeter 
amplitude as ^ pi , P2 (xi, x 2 ) = c lPX ip q (x). Then, the integrals over X and X 1 are performed and the formula (|A1[) 
changes into 



C(P,q) = / d 4 xd*x'pp(x;x')ip q {x)ip* q (x') . 



And now one argues that the density matrix pp(x; x') can be approximated by the diagonal form 

PP {x-x') = D r {x) 8^{x-x') , 



(A2) 



(A3) 



where D r (x) is the relative source function of probabilistic interpretation. The nonrelativistic counterpart of D r (x) is 
given by Eq. ([7|). To justify the expression (|A3j) one assumes that the effect of particle production can be factorized 
from the final state interaction, as the production process occurs at a much larger energy-momentum scale than the 
process of final state interaction. Substituting the formula (|A3|) into Eq. (|A2|) one finds 

C(P,q) = J d A xD r (x) \i> q {x)\ 2 . (A4) 

When the Bethe-Salpeter amplitude 4> q (x) is transformed to the center-of-mass frame of the pair of particles, it can 
be replaced by the nonrelativistic function y q , (r*) when t* is assumed to vanish. Then, one reproduces our formula 



APPENDIX B: APPROXIMATION OF CONFLUENT HYPERGEOMETRIC FUNCTION 



We derive here an approximate expression of the Coulomb scattering function which holds when the source size is 
much smaller than the Bohr radius of the two interacting particles. The confluent hypergeometric function i-Fi(a, b; z), 
which gives the Coulomb scattering function, is defined as 



°° z" Vr 1 a + k 



1 F 1 (a,b;z) = l+ -r II 



m b + k 

n = 1 k = 



(Bl) 



The Coulomb scattering function corresponds to the arguments a — —irj/q, 6=1 and z = i(qr — qr). Introducing 
the parabolic coordinate £ we have q£ = qr — qr and thus, 



y n=l y ' k=0 H 



(B2) 



To obtain the desired approximation we write down a few first terms of the series (|B2j) and we rearrange them as 



/ in \ ~-iqZ - H (- H + l)(«^) 2 



(3!) 2 



17 



1.30 
1.25 
1.20 
1.15 
1.10 
1.05 
1.00 
0.95 
0.90 



exact 

approximate 




-4 = l fm 



10 20 30 40 50 60 70 80 90 100 

q [MeV/c] 



0.16-1 lm[ 1 F 1 ] 
0.14 
0.12^ 
0.10 
0.08 
0.06-\ 
0.04 
0.02 
0.0Q-\ 



exact 

approximate 




10 20 30 40 50 60 70 80 90 100 



q [MeV/c] 



FIG. 15: The real (left panel) and imaginary (right panel) part of the pion Coulomb function iFi^—irj^/q, l;iq£) as a function 
of q for several values of £. The solid and dashed lines represent, respectively, the exact formula (|B2|) and the approximation 

(TBI. 



V)2 



(l!) 2 (2!) 2 (3!) 

oo - n— 1 
n=i y '' fc = 



We first note that with the formula (|B3j) . one easily finds the value of the Coulomb function at q = 

lFl ( - 7' 15 iq t) L = 1 + E W II (»*) =ESf = *>(2>/tf ) , (B4) 

y V n=l V fc = n = y '' 

where I v (z) is the modified Bessel function of the first kind defined as 

^ 1 / z \2k+v 

Iv{z) = E T (k + v + i)k\ (.2) 

We define the new variables x = r/t; and y = iq/rj, and we write down the series (|B3p as 

P /l , \ x , x 2 (l + j/) x 3 (l+y)(l + 2^) 



.y' ' tf / (l!) 2 (2!) 2 (3!) 2 

n— 1 n— 2 n— 3 

We are interested in the approximation which holds when the source size is much smaller than the Bohr radius of the 
scattering particles that is when x <C 1. Since y can be arbitrary big, the series cannot be simply terminated at a 
given power of x. Instead, one should take into account the lowest power of x for every power of y. For this purpose 
we have to rearrange the series (|B5|) . After rather tedious analysis, one shows that 

oo oo / n—k n— fc+1 n— 1 \ n 

^4H = 1 + I> fe E E E - E ot- (B6) 

y fc = n = fc+1 \Zi = lZ 2 =Zi + l Z fc =Zi+fc-l / v ' ; 



And now for each k in the series (|B6|) we take into account only the term of the lowest order of x that is we include 
only the term of n = k + 1. Observing that 

12 fc 

EE- E iife. ••** = *!, 

ii = l/ a = Ji+l i k = i 1 +k-l 
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FIG. 16: The real (left panel) and imaginary (right panel) part of the kaon Coulomb function iF\{—irjK / q, 1; iqQ as a function 
of q for several values of £. The solid and dashed lines represent, respectively, the exact formula (|B2|) and the approximation 
(IB81). 



we obtain the desired approximation 



Fi(-,l;xy)*l+-Y l 



(xy) 



fe+i 



l 



/ , , , , . — rr = 1 + Si(-ixy) (7 e + ln(-iary) - Ci(— ixy)) 



(B7) 



where Si(z) and Ci(z) are integral sine and cosine functions, respectively, and j e f=a 0.5772 is the Euler's constant. 
Reintroducing the physical arguments, we finally have 



if! ( - ^, 1; igC) « 1 + ^Si(gO + * ^ + ln( 9 - C%0) 



(B8) 



In Figs. [TS] and \W\ we show the Coulomb functions computed from the approximate formula (|B8[) for the Bohr 
radius of pions and kaons, respectively. As seen, the approximation works very well for pions {ri~ l = 388 fm) but 
it is not so accurate for kaons (ry^- 1 = 110 fm). We also see that Re[i-Fi(— ir)/q, 1; iqC)] ^ ir]/q, 1; iq(,)]. 

Therefore, the imaginary part can be neglected and 



,F,.( - l l,l;iqt) »l + 2si(g0 
q J q 



Since ^Si(g£) C 1, we also have 



i^i ( 



l;iq£ 



1 + 2^Si(gO 

q 



(B9) 



(B10) 



The approximations (|B9l IB10[) were used to compute the correlation functions of pions but not of kaons. 
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